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In this work, a cosmological model inspired by string/M- theory with fermionic field is 
taken into consideration. Here it is investigated whether the introduction of a non-Dirac 
fermionic field— characterized by an interaction ierm-affects the cosmological evolution. 
The self-interaction potential is considered as a combination of the scalar and pseudo- 
scalar invariants. It is observed that the fermionic field under consideration behaves like 
an inflation field for the early Universe and later on, as a dark energy field. The late time 
acceleration becomes more prominent by the addition of the interaction term. There is 
a slight decrease for the inflation peak as well as for the energy density. We see that 
addition of higher order terms to the fermionic part of Lagrangian does not significantly 
change either the inflation or the late time acceleration behavior. 
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1. Introduction 

In cosmology, the investigation for the constituents responsible for the accelerated 
periods in the evolution of the universe is of great interest. The mysterious dark 
energy has been proposed as a cause for the late time dynamics of the current 
accelerated phase of the universe. Recently several approaches were made to explain 
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accelerated expansion by choosing fermionic fields as the gravitational sources of 



Ribas et.alW, in 2005, investigated whether fermionic field, with a self-interacting 
potential that depends on scalar and pseudo-scalar invariants, could be responsible 
for accelerated periods during the evolution of the Universe, where a matter field 
would answer for the post inflation decelerated period. They have shown that the 
fermionic field behaves like an inflation field for the early Universe and later on, 
as a dark energy field, whereas the matter field was created by an irreversible 
process connected with a non-equilibrium pressure. They observed that for the old 
decelerated phase of the universe the fermionic field plays the role of dark energy 
and drives the universe into an accelerated regime. But it is extremely difficult to 
present a completely satisfactory theory of dark energy. Hence cosmologists now 
look to string/M-theory to formulate a successful dark energy theory because string 
theory could provide some insight on the peculiar initial conditions required by 
standard big bang. Nojiri et.al ™ proposed the Gauss-Bonnet dark energy model, 
inspired by string/M-theory where standard gravity with scalar contains additional 
scalar-dependent coupling with Gauss-Bonnet invariant. Their study indicates that 
current acceleration may be significantly influenced by string effects. 

In the present work, we propose a cosmological model with fermionic field. Notic- 
ing that string/M-theory requires additional higher order terms, we introduce pow- 
ers other than those that would emerge from the potential of Ribas et.al. To do so 
in a systematic way, a term proportional to the square of the fermionic Lagrangian 
density, which we call as an interaction term, is added to the total Lagrangian. 
It may be noted that the fermion condensate expected from super symmetry also 
calls for higher order terms. In the present paper our aim is to investigate how the 
non-linear interaction term affects the cosmological evolution. 

2. String Theory Inspired Cosmological Model with Fermionic 
Field 

The action for this model reads 



where L g = R/2, with R denoting the curvature scalar, L g is the gravitational 
Einstein Lagrangian density, L m is the Lagrangian of the matter field and di and 
ci2 are coupling constants. We are working under the assumption that higher powers 
of Lq grow less important. The L 2 D term is suppressed relative to the main terms 
by the fourth power of the Planck mass, M v i . (The Lfj term would have eighth 
power of M p i suppressing it relative to the main term etc). The values of coupling 
constants a\ and a-i cannot be bigger than unity for the same reason; otherwise the 
assumption of ordering would be wrong. We have M~ t 2 = 8itG . Later, this is set to 
be equal to unity. Finally, the Dirac Lagrangian density Ljj of the fermionic field 



energy. 




(1) 



December 7, 2009 15:49 WSPC/INSTRUCTION FILE Rakhi 



A Cosmological Model with Fermionic Field 3 



HOI 



for a fermionic mass to is given by 



In Eq. |(2]), = V(\l/, ^) describes the potential density of self-interaction between 
fermions and W = ^^7° denotes the adjoint spinor field. Moreover, the connection 
between general relativity and Dirac equation is done via the tetrad formalism and 
the components of the tetrad play the role of gravitational degrees of freedom. That 
is, T' 1 = ef^y a are the generalized Dirac-Pauli matrices, where denote the tetrad 
or vierbeiri^. The covariant derivatives in Eq. |2]) are given by, 

D^ = d^+m^ (3) 



with the spin connection, 



1 



% = --g pa [T^-e p b d,e b 5 ]T°r* (4) 
where T^ s is the Christoffel symbol. 

2.1. Basic Equations 

Through Euler-Lagrange equations, from Eqs. and |2j), we can obtain the equa- 
tions of motion for the spinor field as 

^ ai +2^L D JdL D /dlp = (5) 

which requires that either ^a,i + = or dLo/dip = 0. 

The first case, i.e. (a\ + — , gives Lp = — -, — a \ r . This appears 

V M pl J 2 { a 2/ M pl) 

to be an unphysical constraint, and is not explored here further. 
We take dLu /dip — as our ansatz. I.e. 

iT^D^ - to* - = (6) 

The equation of motion is exactly the same as in the case where no L 2 D term is 
present. However we see that the L 2 D term modifies the source terms for gravity 
in a significant way. We also note that any series of arbitrary powers LJ, present 
in the original Lagrangian would lead to SL^/Sip = for the fermionic equation 
of motion, because the other equation arising would be unphysical. So, with the 
above ansatz, our analysis can be thought of as treating the case of a general form 
for the fermionic contribution to the stress tensor, where we retain terms only upto 
the next to leading order contribution in an expansion of the stress tensor in the 
derivatives of i/>. 
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Similarly, we can obtain the Dirac equations for the adjoint coupled to the 
gravitational field. 

The variation of the action Eq.([T]) with respect to the tetrad gives the Einstein's 
field equations as 

R/j,v - -^g^uR = —T^v (7) 

where — + , with Tj being the energy- momentum tensor of the 
fermionic field and Tfif, that of the matter field and T M „ = — -^=^|j- . The sym- 
metric form of the energy-momentum tensor of the fermionic field which follows 
from Eq. ^ gives 

f 4 

-L D {l>T fl D v iS + ^VD^ - D U 9T^^ - D^T u ty}- 

g^ ai Ln-g^^Ll (8) 
m pi 

where the Dirac matrices have the form F° = 7 , T l = ^75-7*, T 5 = 7 5 and the spin 
connection components are 

n = 0,O i = ia(t)7 i 7° (9) 

Based upon the cosmological principle, which states that our universe is homogenous 
and isotropic, we use the Robertson- Walker metric to describe our Universe: 

ds 2 = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 ) (10) 

The interactions between the constituents are modeled through the presence of a 
non-equilibrium pressure term (vu) in the source's energy-momentum tensor. There- 
fore, for an isotropic, homogenous Universe, the total energy-momentum tensor- 
which is composed by fermionic and matter fields- can be written as 

T" = diag(p, —p — vj, —p — w, —p — vo) (11) 

where p = pf + p m and p = pj +p m and the quantity w refers to a non-equilibrium 
pressure which is related to dissipative processes during the evolution of the universe 
and represents an irreversible process of energy transfer between the matter and the 
gravitational field. By the use of Eq. (|lip . the components of the energy-momentum 
tensor of the fermionic can be identified as 

= diag{p f: -p f ,-p f , -p f ) (12) 

Now the non-vanishing components of the energy-momentum tensor for the 
fermionic field following from Eqs. fl}, and (fTQ|) can be written as, 

PWS = * (m <«> + V)-^ (m {«> + V)\^ (™ <«> + |f + f f 

(13) 
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(Tf)\ = (T f )l = (Tf) 



a i 



V 



fdV 
~2 df 



d¥"2 



a-2 



V 



^dV 
~2 df 



d¥"2 



Now, the self-interaction potential takes the form, 

\ 2 



V = A 



(14) 



(15) 



where A is the coupling constant and n is a constant exponent. We consider V 
as a combination of the scalar and pseudo-scalar invariants. We assume both the 
scalar and pseudo-scalar densities to play equally important roles and assume the 
coefficients j3% and to be both O (1). For simplicity we consider the values to be 
unity. That is, j3x = 02 = 1. 

Combining Eqs. (fT2")) . (|T5|) . (TTJ| and (TT5|) we solve for the energy density and 
pressure of the fermionic field as 

Pf = ai (to (m) + V) - ^ (m (tftf ) + + (m + 2n^) 2 (16) 



CI 2 



P/ = Oi (2n-l)V+-^[(2n-l)V 



A/f 4 



(17) 



respectively. We infer from Eq.([T7|) that the fcrmions could be classified according 
to the value of the exponent n. Indeed, for n > 1/2 , the fermions represent a matter 
field with positive pressure (n > 1/2) or a pressure-less fluid (n = 1/2), whereas for 
n < 1/2 the pressure of the fcrmions is negative and they could represent either the 
inflation or the dark energy. For the massive case, we shall deal with only the case 
where the fermionic field behaves as inflation or dark energy, i.e., the case where 
n < 1/2.12 

The Friedman equations are 

9 8ttG „ n 8nG 8ttG , 
ZH 2 = — T ° = — p = —{p f + p m ) (18) 



2 a - = 
a 



8ttG 
Ik 2 " 



T ° + 3Tl] 



or 



a 



BttG 



[Pf + Prn + 3p f + 3p m + 3w] 



a 3c 2 

Using Eqs. (JTHJ) and (|2ll|) . the acceleration equation becomes 

a _ 8ttG 
a 6c 2 



[Pf 



3p f + 3p m + 3w] 



(19) 



(20) 



(21) 



The covariant differentiation of Einstein field equations, using the Bianchi identities, 
leads to the conservation law of the total energy-momentum tensor: i.e T£ v = 0. 
I.e. 



P+3H(p + p + m) = 



(22) 
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where H = 2^ is the Hubble's parameter. The conservation law for the energy 
density of the fermionic field and matter field can thus be written as 

Pf+3H(j> f + Pf )=0 (23) 

p m +3H (p m + p m + w) = (24) 

Considering the fermionic field to be a function of time alone, the Dirac equations 
Eq.© becomes 

q>+-H^ + im-f°y + ij ^ = (25) 
2 d<J/ 

Similarly, the Dirac equation for the adjoint spinor field coupled to the gravitational 

field becomes 

~$+-H~$-im^° -i^7° = (26) 

2.2. Field Equations 

The system of field equations that help us to find the cosmological solutions is: 



2.2.1. The Acceleration Equation 

The acceleration equation for the cosmological model under consideration can be 
obtained from equations (|16[) . (JT7J) and (f2"Tj) .We have chosen the units so that h = 
c = 1. Also we consider the pressure of the matter field as p m — w m p m with 
< w m < 1, which is a barotropic equation of state. We again consider a vanishing 
energy density of the matter field (p m ) and a vanishing non-equilibrium pressure 
(n7)at t=0. Incorporating all these ideas, the normalized acceleration equation for 
the present model becomes 

a = _8tt Pm + 3 pf + 3p m + 3w) (27) 

a 6 

where 

Pf = at (to (!^) +V) - jft- (m (!^) + V) 2 + -$r (m <*^) + 2nV) 2 
p f = ai (2n - 1) V + a 2 [(2n - 1) V] 2 



Putting 8ttG = 1 , Eq. ||27J will become 
d 1 

- = -g (Pf + Pm + 3p f + Zpm + 3tu) (28) 

with 

p f = ai (to + V) - a 2 (m + V) 2 + a 2 (m + 2nV") 2 

P/ = a t (2n -l)V + a 2 [(2n - 1) F] 2 



For further calculations we use this result. 
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2.2.2. Evolution equation for the energy density of matter field 
The evolution equation for p m is given by 

p m +3H(p m + W m p rn + w) = 



(29) 



2.2.3. Evolution equation for the non- equilibrium pressure 

The linearised form of the evolution equation for non-equilibrium pressur j^EEl reads 

r w +m + 3?]H = (30) 

where the coefficient of bulk viscosity 77 and the characteristic time t are assumed 
to be related with the energy density p by 77 = ap and r = - , where a is a constant. 



2.2.4. The Dirac Equation 

In terms of the spinor components,^ = ^2, ^3, ^ / 4) T , the Dirac Eq. (|25|) can 
be written as 



d 
di 



#2 
*3 

v* 4 / 



#2 



-2i + *1*2 - *1*3 - *2*4) 



/*3 \ 



d(*7 5 *)' 



= 



*2 
-*3 

(31) 



2.3. Cosmological Solutions 

In order to obtain numerical solutions of the coupled system of Eqs. ([27| - (|3"Tj) . we 
have to specify the parameters A, /3i, /?2, w, a, w m and a 1 . These parameters take 
the following values: 

A = 0.1, j3\ = fa = 1, n = 0.25, m = 0.01, a = 1, u> m = 1/3 and aj = 1. 



2.3.1. Dependence of the evolution on the coefficient of the interaction term 
a 2 

The initial conditions, we have chosen for t=0 (by adjusting clocks) are: 
a (0) = 1, tf 1 (0) = O.li, 4-2 (0) = 1, * 3 (0) = 0.3, 4< 4 (0) = i, p m (0) = 0.05, w (0) = 
0. The value of p m (0) is so chosen that the observable normalized matter density is 
5%. The last initial condition corresponds to a vanishing non-equilibrium pressure 
at t = 0. The conditions chosen here characterize qualitatively an initial proportion 
between the constituents in the corresponding era. 
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Graphs are plotted for various values of et2. The effect of L 2 D is now studied first 
by choosing a 2 = 0.5 [Figs. 1-2]. 

In Fig. 1, we plot the evolution of accleration versus time for a 2 = 0.5 while Fig. 3 
shows the the evolution of accleration versus time for a 2 = 0. In Fig. 2, the evolution 
of energy densities of fermionic field and matter field is plotted. 
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Fig. 1. Acceleration vs. time ai = 0.5 




Fig. 2. Energy density of fermionic pf and matter p m fields vs. time 




Fig. 3. Acceleration vs. time for ci2 = 0[Ribas et.al(2005)] 
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The effect of L 2 D for various values of 02 is studied below by taking 02 = 
0.1,0.2,0.8 respectively [Figs. 4-5]. 




Fig. 4. Acceleration vs. time for various 012 
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Fig. 5. Energy density of fermionic pj and matter p m fields vs. time for various 

In Fig. (1) and Fig. (4), acceleration vs. time is plotted whereas in Fig. (2) 
and Fig. (5), it is shown the behavior of the energy densities of the fermionic pf 
and matter p m fields as a function of time . We infer from these figures that the 
proposed model behaves like an inflation field for the early Universe and later on, 
as a dark energy field. Here, the late time acceleration becomes more prominent 
by the addition of the interaction term. There is a slight decrease for the inflation 
peak. 

2.3.2. Study of evolution for (i)tp/M^A « 1 and ($ip/M^ » 1 

We have already mentioned that the power n of the scalar and pseudo-scalar den- 
sities in EqfT5]is required to be less than 0.5 for the fermionic field to behave like a 
dark energy at late times. With the addition of the interaction term, n will have to 
be less than 0.25 to avoid suppression of this term. Then if Otp'tp / << 1 (i.e. 
< 1 in Eql28|. the relative contribution from the V in L|,will be larger compared 



December 7, 2009 15:49 WSPC/INSTRUCTION FILE Rakhi 



10 Rakhi et al. 

to the case when ^ip > 1. Thus the new effects we obtain with the addition of the 
interaction term should be more prominent when scalar or pseudo-scalar density is 
less than one, and when the densities become greater than one, the results obtained 
with the addition of the interaction term should be closer to the behavior seen with 
Ld alone. 

In order to check this, we plot the acceleration versus time, when scalar or 
pseudo-scalar density (at t=0) is greater than one [Fig. 6], and less than one [Fig. 7] 
, for a 2 = and for a 2 = 0.5 in each case. 

The initial conditions, we have chosen for tptp < 1 at t =0 are: a (0) = 1, "fi (0) = 
O.li, * 2 (0) = 1, * 3 (0) = 0.3, * 4 (0) = i, p m (0) = 0.05, w (0) = 0. 
whereas the initial conditions chosen for 'ijjip > 1 at t =0 arc: a (0) = 1, (0) = 
O.Oli, * 2 (0) = O.Oli, * 3 (0) = 1, * 4 (0) = i, p m (0) = 0.05, w (0) = 0. 
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Fig. 6. Acceleration Vs time for ipip < 1 
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Fig. 7. Acceleration Vs time for tptp > 1 



From the figures, we notice that the change in the accelerating behavior produced 
by the addition of the interaction term is more prominent, when scalar or pseudo- 
scalar density (at t=0) is less than one. In comparison we note that the difference 
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is quite small when scalar or pseudo-scalar density (at t=0) is greater than one, as 
expected. 

2.4. Final Remarks and Conclusions 

In the present work, we have considered a cosmological model inspired by string/M- 
theory with fermionic field. Here we have assumed that there exists a field with a 
kinetic term like L 2 D which we call as an interaction term and this field corresponds 
to dark energy which contributes to expansion/inflation etc. The self-interaction 
potential is considered as a combination of the scalar and pseudo-scalar invariants. 

We have investigated the contribution to the cosmological evolution of a non- 
Dirac fermionic field-characterized by an interaction term L 2 D . Early universe the- 
ories require the presence of higher order terms. To study these in a systematic 
way, we have included the L 2 D term in the total Lagrangian and checked the effects 
of this non-linearity. We see that the model behaves like an inflation field for the 
early Universe and later on, as a dark energy field, similar to what Ribas et.al. 
obtained [Fig. 3] with the addition of Lr> only. Our results show that the higher or- 
der non-linear terms do not significantly change the overall behavior, obtained with 
Ld alone; except that the late time acceleration is even more prominent. Also it is 
significant that the change in the accelerating behavior produced by the addition 
of the interaction term is more prominent, when scalar or pseudo-scalar density (at 
t=0) is less than one, whereas the difference is quite small when scalar or pseudo- 
scalar density (at t=0) is greater than one. So we conclude that fermionic terms 
to all orders in the Lagrangian lead to inflation at early times and acceleration at 
late times for the universe. Our model can generate standard acceleration similar to 
Dirac behavior as well as a non-Dirac behavior by simply adjusting the magnitude 
of the coupling constant a2 . 
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